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ABSTRACT 



A method using electrical resistivity measurements to deter- 
mine the in situ porosity of marine sediments was investigated in 
the laboratory, and equipment for this purpose was designed, fabri- 
cated and tested. Formation Factor-porosity relationships deter- 
mined in the laboratory for three clays (kaolinite, illite, and 
montmorilloni te) , Providence silt, four sands, and four marine sedi- 
ments showed that porosity was predicted within + 2 percent. The 
Formation Factors ranged from 1,1 to 5,9 while porosity ranged from 
26 to 93 percent. The particle size and distribution influenced the 
electrical resistivity of these sediments independent of porosity 
while particle shape did not. The laboratory equipment was econo- 
mical, safe, easy to operate and could be used to determine perme- 
ability, tortuosity and void ratio-log consolidation pressure of 
sediments • 

The Electrical Resistivity Measuring System designed for the 
Deep Ocean Sediment Probe consists of three interchangeable electrode 
arrays, the electronic circuit and the FM telemetry data link with 
6000 feet of coaxial cable. The predicted porosity values obtained 
with the System in the laboratory were within + 2 percent of the line 
of best fit obtained with the earlier laboratory equipment. Using 
the inner corer ring electrode array, the System was tested in the 
shallow water of Narragansett Bay and most of the data agreed well 



with the value obtained later with the laboratory equipment on the 
same sediment* 

The inner corer ring electrode array read continuous conducti- 
vity through the water/sediment column thus permitting a value of 
salinity to be obtained* The bottom water salinity obtained in this 
way was in good agreement with the interstitial water salinity of the 
cores. The system detected a core loss, touch down, and sample dis- 
turbance. This relatively low cost system allows accurate, quick 
in situ porosity determination and shows promise in other areas such 
as the monitoring of pollution, the prediction of other bottom sedi- 
ment properties and as a warning system on ocean structures. 
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PREFACE 



This dissertation is organized according to the standard plan. 
The text is divided into two sections. The first section includes a 
report of laboratory experiments on an indirect method of determining 
porosity of sediments and a review of the literature involving similar 
work. The instrumentation used in the laboratory is described in 
this section. 

The second section describes the design, fabrication, testing, 
and evaluation of an in situ system used to determine porosity of 
marine sediments and the equipment used by others for similar work. 
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I 

THE USE OF ELECTRICAL RESISTIVITY TO PREDICT 



THE POROSITY OF SATURATED SEDIMENTS 



INTRODUCTION 



The purpose of this work is to investigate the accuracy of 
using an electrical resistivity method to predict the porosity of 
soil-water mixtures. Porosity in a two phase soil water mixture is 
the ratio between the volume of water to the total, volume. The deter- 
mination of porosity is important to ocean engineering projects such 
as underwater acoustics, foundation engineering, submarine soil me- 
chanics and bottom survey work. At present, porosity is determined 
in the laboratory. If electrical resistivity is found to be an accu- 
rate method of predicting porosity, a practical method of in situ 
testing would be available. This method would be quicker, more eco- 
nomical, and allow more data to be taken than the present method. 

To verify the accuracy of the electrical resistivity method, 
extensive laboratory tests we re conducted and the results were com- 
pared with theoretical and empirical work reported in the literature. 
Most of the earlier work studying the relationship between electrical 
resistivity and porosity was conducted on formation rocks or disper- 
sions. The geophysical techniques used in well logging analysis 
involve electrical resistivity measurements in order to obtain: 

1. Subsurface correlation with formations in other wells, 
structural mapping. 

2. Lithology of the formations. 

3. Depth and thickness of productive zones. 
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4 . Formation porosities. 

5. Oil and gas saturation. 

6. Volume of oil in a reservoir. 

The porosity of formation rocks ranges from 5 to 30 percent, 
while the porosity of emulsions and dispersions ranges from 50 to 
90 percent. Although these ranges cover the porosity values of 
sediments (which range from 26 to 90 percent), extensive testing is 
necessary before it can be considered valid to apply the models and 
theories used for rocks, emulsions, and dispersions to sediments. 

Electrical resistivity has also been used in order to deter- 
mine the dielectric constant of pure crystals, the volume concentra- 
tion of certain biological materials such as red blood cells, and 

1 2 

ground water detection. (Meredith , Patten and Bennett ) 

Very little electrical resistivity work has been conducted on 
sediments. The limited work that has been conducted on sediments 
generally has not considered clay mineral content and type, particle 
size, shape, distribution and structure, and consolidation pressure. 

In this study the following parameters were varied in order to 
measure their influence on the prediction of porosity as determined 
by electrical resistivity measurements: 

1. Type of sediment (sands, silt, kaolinite, illite, and 
montmorillonite clays, and marine sediments), 

2. Size of particles, 

3. Shape of particles, 

4. Distribution of particle size, 

5. Structure of particles, 
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6. Consolidation pressure, 

7. Salinity of interstitial water. 

The temperature was measured and a correction factor was applied in 
order to standardize all electrical resistivity measurements to 
25° Centigrade. By measuring the electrical resistivity as each 
parameter was varied, while at the same time determining porosity by 
the standard method, the influence of each parameter on the accuracy 
of this method of porosity prediction could be evaluated. 

Since the terms found in this work are used both in submarine 
soil mechanics and geophysics, they are defined below in order to 
eliminate confusion. 

The two-phase mixture of soil and water is considered a satu- 
rated sediment, or simply a sediment. The electrical resistivity 
measurements are taken on both the sediment and the interstitial 
water. The resistivity of the sediment is divided by the resistivity 
of the interstitial water to obtain a ratio called the Formation 
Factor. This ratio is used as a direct index of porosity. 

Electrical resistivity (specific resistance) is a measure of 
the ability of a solution, solid or mixture of unit length and cross 
section to resist the flow of electric current. In solids, current 
is carried by free electrons; in solutions, the current is carried by 
ions. As electrons and ions move through solids and solutions, a 
frictional drag results. This frictional drag per unit volume is 
considered the electrical resistivity of the solid or solution. 
Electrical conductivity (specific conductance) is a measure of the 
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ability of a solution, solid, or mixture, of unit length and cross 
section to carry an electrical current. Consequently resistivity 
and conductivity are inversely related measures of the same property. 
Studies dealing with formation rocks have related resistivity to 
porosity while studies dealing with emulsions and dispersions have 
usually related conductivity to porosity. 

Listed below are some of the symbols, assumptions, and defini- 
tions used in this paper: 

dispersed phase (d) - solid soil particle, 
continuous phase (c) - interstitial water, 
mixture (m) - the mixture of the continuous and dispersed 
phases - sediment. 

The properties of these phases are assumed to be isotropic and uni- 
form. 

r , - resistivity of the continuous phase dispersed 

c , q , m 

phase and mixture respectively, 

k , - conductivity of the continuous phases dispersed 

c , u , m 

phase and mixture respectively. 

In order to make relationships dimensionless the resistivity 
and the conductivity of the continuous phase is considered unity. 
These dimensionless notations are denoted by Roman capital letters: 



r 




t 



( 2 ) 
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_ m 

R = — , 
m r f 
c 



(3) 



k = = i , 

c k f 



(4) 



K d = r ■ 

c 



(5) 



k 

K = 
m k 

c 



( 6 ) 



Formation Factor (FF) which is often used in the Geophysical 

literature is equal to R and the reciprocal of K : 

1 m r m 



FF = R = ~ 
m K 



m 



(7) 



Porosity is normally defined as the volume of the voids to the total 
volume. Porosity in a saturated sediment (two phase system) can be 
defined as the ratio of the volume of water to the total volume: 



_ * a. / \ voids 

Porosity (n) =* - f 

total 



( 8 ) 



if saturated 



V . 
water 

^total 



<4 

( 9 ) 



REVIEW OF THE LITERATURE 



The literature contains many theoretical and empirical models 
that were developed in order to predict porosity from electrical 
measurement. The following review examines the models that might 
logically be applied to a saturated sediment. No one model can 
cover the wide range of sediment types and mixtures (sands, silts 
and clays). Factors such as particle size, shape, distribution and 
structure, and consolidation pressure influence model selection. 

In addition conduction of the clay fraction changes as the ionic 
concentration of the interstitial water varies. The first part of 
this review lists some of the properties peculiar to sediments that 
need to be considered in models in order to enable accurate porosity 
prediction from electrical measurement. 

Except for native metals, sulfides, a few oxides, graphite, 
and high-grade coal, the resistivity of most minerals that make up 
sediments is measured in millions of ohm-meters. Although the soil 
mineral is non-conductive, the soil particle carries an electrical 
charge. The net electrical charge may arise from any one or a combi- 
nation of the following factors: 

1. Isomorphous substitution (important in illite and 
montmorilloni te clays), 

2. Broken bond charges (important in kaolinite clays), 

3. Exposed hydroxyl groups, 
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4. Absence of cations in the crystal lattice, 

5. Presence of organic matter* 

The cation exchange capacity is a measure of the electrical 
charge carried by a soil particle. The magnitude of the cation ex- 
change capacity of a sediment depends on the amount and type of clay 
fraction present and the clay particle surface area® The clay frac- 
tion is usually defined as that part of the soil mass which has an 
equivalent Stokes diameter of less than two microns. Particles 
smaller than one micron are considered in the colloidal range. It 
is in the colloidal range that surface properties begin to dominate 
the physical and chemical properties of materials. 

For most clay particles the specific surface (the magnitude 
of the surface area per unit mass) is a good indicator of the rela- 
tive influence of electric forces on the behavior of the particle. 
Table 1 shows the cation exchange capacity and specific surfaces of 
the common clay minerals. 
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TABLE 1 



THE CATION EXCHANGE CAPACITY AND SPECIFIC SURFACE 
OF KAOLIN ITE, ILLITE, AND NONTMOR ILLONITE CLAYS 



CLAY MINERAL 


^100 g ; 


2 

<7-> 

g 


CATION EXCHANGE CAPACITY ' 


*v> 

SPECIFIC SURFACE 


Kaolini te 


3 - 15 


10 


Illite 


10 - AO 


100 


Montmorillonite 


80 - 150 


1000 



VC VC Cf 

Lambe and Whitman 

When individual clay particles are dropped into water the ions on 
the surface of the particle tend to move out from the surface to be- 
come exchangeable ions. Water molecules are attracted to these ions 
and the two together form what is called the double layer around the 
clay particle. The thickness of the double layer is defined as the 
distance from the soil surface to the point where the ion concentra- 
tion is in equilibrium with the interstitial water. 

Many factors influence the size of the double layer. The 
Gouy-Chapman Double Layer theory is often used to explain the double 
layer and the factors affecting it. This theory assumes that there 
is no interaction between clay particles and that the clay particles 
are plates. This theory states that the size of the double layer is 
directly proportional to the temperature and dielectric constant of 
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the solution and inversely proportional to the concentration and 

ionic valence of the electrolyte in the interstitial water and the 

* 

electric charge on the clay particle. Other factors influencing 
the double layer are pH and size of ions in solution. 

If the cation exchange capacity and the specific surface are 
low and a small double layer exists, almost all the interstitial 
water would be free water and would not be attracted to the mineral 
skeleton. Conversely if the cation exchange capacity and the specif- 
ic surface are high and a large double layer exists, it would be 
quite difficult to separate the clay mineral and interstitial water 
phases . 

One of the earliest works examining the effect of charged 
clay particles and the associated double layer on electrical resis- 
tivity was done on reservoir rocks by Winsauer et al.^ These in- 
vestigators concluded that charged clay particles do influence the 
conductivity of shaly sand and limestone formations, but that this 
influence became negligible at high ionic concentrations of inter- 
stitial water. 

After conducting laboratory experiments, Hill and Milburn^ 
made similar conclusions and developed empirical equations for 
reservoir rocks relating Formation Factor to the resistivity of the 
interstitial water, porosity and the cation exchange capacity. In 
their equations as the electrical resistivity of the interstitial 
water decreased, i.e., became more saline, the effect of the cation 
exchange capacity decreased. 

One of the first works to study the effect of the cation ex- 
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change capacity on the electrical resistivity of a sediment was con- 
ducted by Berg. ^ Using kaolinite clay and sodium chloride solutions, 
he showed that the Formation Factor remained constant when the water 
resistivity was less than 80 ohm- cm. However when water resistivity 
was greater than 80 ohm-cm, the Formation Factor decreased. Water 
that has an electrical resistivity of 80 ohm-cm would have a salin- 
ity of seven parts per thousand. Sea water having a salinity of 

35 parts per thousand would have a resistivity of 25 ohm-cm. This 

5 6 

work agrees with the work on reservoir rocks. f The effect of the 
charge of kaolinite clay particles on the electrical resistivity ap- 
peared to be negligible at high ionic concentrations of interstitial 
water. 

g 

Sheeler et al. expanded the knowledge on sediments by study- 
ing several types of clay. He found that for kaolinite clay the ef- 
fect of cation exchange capacity on the electrical conductivity and 
thus Formation Factor persisted up to a 3 percent sodium chloride 
solution, then became negligible. However when using montmor i 1 loni te 
the ion effect persisted up to an 8 percent sodium chloride solution. 
Sea water is approximately 3.5 percent sodium chloride. This study 
suggests that the type of clay present in a sediment such as 
inontmorillonite may influence the electrical resistivity measurements 
due to the cation exchange capacity even though the salinity is 

higher than sea water. 

9 

Boyce reported that the percentage of clay size material 
varied inversely with the electrical resistivity of fine grained 
marine sediments obtained from the Bering Sea. Boyce explained these 
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results could be due to: 

1. a specific conductivity increase resulting from a 
double ionic layer of clay mineral in the clay size 
fraction, 

2. a secondary conductivity increase caused by an increase 
in porosity, since the clay size fraction and porosity 
are directly related, 

Fritsch and Tauber^ studied kaolinite, illite, and 
montmorilloni te clays at low and high salt concentrations. They 
observed that the electrical resistivity measurements were over- 
whelmingly influenced by the specific surface area of the clay par- 
ticles when distilled water was used, but at high salt concentrations 
the resistivity measurements were independent of this influence. As 
illustrated in Table 1, the specific surface area is an indicator 
of cation exchange capacity. These authors concluded that further 
work was necessary in order to establish the limits of this influ- 
ence for different clays at different concentrations, 

Atkins and Smith*’*' studied sodium montmorilloni te, calcium 
montmorilloni te and kaolinite, which have cation exchange capacities 
of 74,0, 60.1 and 16.1 meq/100 g respectively. They found that the 
Formation Factor remained constant even though the interstitial 
water conductivity ranges x^ere 1.1 to 8.7, 0.1 to 19.7, and 0.009 
to 8.5 mho/meter respectively. They found the same result x*ith 
illite clay, but the range of interstitial water conductivity used 
was not given. 

12 

Kermabon et al. further studied the effect of ionic solu- 
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tions on the cation exchange capacity of clay particles. In their 
laboratory experiments the marine clay particles had a grain size 
smaller than 3,9 microns and the sodium chloride solutions ranged 
from 0,7 to 10 percent. The porosity of the clay sediment samples 
was constant. Throughout this wide range of water salinity the re- 
lation between the conductivity of sediment and the conductivity of 
the interstitial water remained constant; therefore, any influence 
the charge of the clay particles had on conductivity was masked by 

the ionic solutions. 

13 

Smith found similar results in laboratory experiments with 
"coarse” clay. He used sodium chloride solutions ranging from 
0.5 to 10 percent. 

The results of these studies indicate that in fresh water the 
charge on clay particles (cation exchange capacity) greatly influ- 
ences the conductive properties of the sediment. This charge and 
consequently its effect on electrical resistivity measurements in 
fresh water depends on the amount and type of clay present. It is 
apparent that in a fresh water environment cation exchange capacity 
must be accounted for by the model or by empirically measuring and 
correcting for its effect on resistivity measurements if accurate 
porosity predictions are to be made. However in a marine environ- 
ment, the influence of the cation exchange capacity of clays on 
electrical resistivity measurements appears to be negligible in most 
cases. The reason for this difference between fresh and salt water 
can be deduced by visualizing the current flow in solution as the 
migration of free ions. In a sediment of clay and fresh water, most 
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of the ions in solution are under the influence of clay particles 
and therefore the current flow through the sediment is inhibited. 

If the interstitial water is saline, there are sufficient ions to 
neutralize the clay particles and allow the current to flow freely. 
If the clay's cation exchange capacity is high, more ions will be 
required to neutralize its effect on conductivity. Although cation 
exchange capacity is usually neglected in the marine environment, 
more testing is needed before this influence on electrical resisti- 
vity measurements for all types of clays can be assumed to be 
negligible. 

The above studies indicate that the cation exchange capacity 
of some clays need not be considered as an influence on electrical 
resistivity measurements in the salinity of ocean water. However 
other characteristics of clay must be considered if accurate poros- 
ity predictions are to be made from electrical resistivity measure- 
ments, The clay type, sand and silt present in a sediment determine 
the size, shape, and distribution of particles, and influence parti- 
cle structure. 

4 

The M.I.T. soil classification of sediment particle size 
is given below: 

from 0,06 mm to 2,0 mm is considered sand, 

0.002 mm to 0.06 mm is considered silt, and 
a particle smaller than 0.002 is considered clay sized. 

The particle distribution of a sediment is defined as a measure of 
spread such as standard deviation from the mean. If the results of 
past work are applied to sediments, it would appear that the parti- 
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cle size and distribution might influence electrical resistivity 
in the following ways: 

If the porosities, particle shapes, and standard deviations 

of two distributions are the same for samples with different mean 

particle sizes, the electrical resistivity will usually be lower 

for the larger particle size. This effect was suggested in cxperi- 

1 14 

ments with glass beads by Meredith and by Wyllie and Gregory. 

As the particle size increased, the electrical resistivity de- 
creased. 

If the porosity and mean particle shape and size are the same 
while the distributions of two sediments are different, the electri- 
cal resistivity will be higher in the'sediment with the greater as- 
sortment of particles. This effect was suggested by studies on rocks 
by Semenov^ and work conducted on glass beads by Meredith.^ 

The effect of shapes in a sediment is complicated. A sedi- 
ment may be made up of particles of different shapes ranging from 
spheres to plates. Usually particles in the sand and silt size range 
are closer to a sphere in shape than the plate-like clay particles. 
The different particle shapes affect the path lengths that migrating 
ions take when current flows. If two sediments have approximately 
the same particle size and spread and equal porosities but differ in 
particle shape, the electrical resistivity would be expected to be 
higher for a sediment composed of disc-shaped particles than a sedi- 
ment made up of spherical shaped particles. This effect was sug- 
gested by studies on various shapes by Wyllie and Gregory .^ 4 Discs, 
cubes, cylinders, and triangular prisms resulted in higher electri- 
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cal resistivity values than spheres at the same porosity values. 

They also found that Ottawa sand (sphericity - 0,895) and glass 
beads (sphericity = 1,00) had similar Formation Factor versus poros- 
ity re lat ionships • 

Meredith^ found that the conductivity of dispersions at a 
porosity of 50 percent resulted in a 10 to 20 percent error if the 
shape difference between spheres and discs was not taken into con- 
sideration. The error decreased with increasing porosity. 

Semenov^ theoretically developed equations for formation 
rocks based on shape. Semenov* s equations for rocks of spherical 
shaped grains and disc shaped grains (having the same distribution) 
resulted in Formation Factor values of 2.7 and 10.1 respectively at 
a porosity of 50 percent. 

Atkins and Smith^ stated that at the same porosity angular 
particles have higher Formation Factors than spherical shaped par- 
ticles. A shape factor was determined for five clays and one mica. 
Shape was inferred to be the major factor causing differences be- 
tween the clays. No electron microscope photographs were taken, and 
data such as mean particle size and distribution were not given. 

In this paper, particle structure will refer to the orienta- 
tion of individual soil particles in the sediment mass. Particle 
structure in sands and silts usually depends on packing. The type 
of packing is important because it affects the porosity and the path 
length required for migrating ions and consequently the electrical 
resistivity measurement. The two extreme types of packing of uni- 
formly sized spherical particles are loose or cubic packing and 
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dense or rhombohedral packing as shown in Figure 1 



FIGURE 1 



SPHERICAL PARTICLES IN A DENSELY AND A 
LOOSELY PACKED CONFIGURATION 





CUBIC PACKING 
(loose ) 



RHOMBOHEDRAL PACKING 
(dense ) 



Structure in clays is greatly influenced by those factors af- 
fecting the size of the double layer. In clays the net electrical 
force between particles is the result of many repulsive and attrac- 
tive forces. If the double layer of clay particles is small, the at- 
tractive forces are dominant and the clay particles will move toward 
each other and become attached (Flocculant Structure). If the double 
layer is large, the repulsive force is more dominant and the clay 
particles will tend to move away from each other (Dispersed Struc- 
ture). The two clay structures are shown in Figure 2, 



FIGURE 2 



ORIENTATION OF CLAY PARTICLES IN A 
FLOCCULANT AND A DISPERSED STRUCTURE 




FLOCCULANT 

STRUCTURE 



DISPERSED 

STRUCTURE 
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In the flocculant structure, large voids may be enclosed by 
the clay skeleton and result in sediments of very high porosities. 

On the other hand the dispersed structure results in lower porosity 
values. Aside from the facts of porosity differences, it would ap- 
pear by visualizing the two structures that the flow of current 
through these two structures would also be different. 

If the current path becomes more tortuous its length is 
increased. Tortuosity of the current path may be defined as the 
ratio of the average statistical length of current flow lines to the 
simple physical length in a straight line from one end of a sediment 
sample to the other. 

By definition the resistance (R) between flat faces of a 
cylindrical sediment sample is: 



R = 





( 10 ) 



where: 

1 is the length of the sediment sample, 

A is the cross-sectional area of the sediment sample, 

1^ is the effective average statistical length of the sedi- 
ment through the pore structure, 

A is the effective average statistical cross-sectional area 
e 

of the pore structure. 



In terms of Formation Factor: 
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1 

e 




c — — 
A A 



where t is the tortuosity factor previously defined. Porosity can 
be expressed as follows:: 



n 



voids 

V 

total 



A 1 



e e 
A1 



A 

e 

A 



t. 



( 12 ) 



If the porosity is the same equations (11) and (12) can be combined 
and tortuosity may be defined as a function of porosity and Formation 
Factor: 



t = */FF*n (13) 

or 

t 2 

FF = — . (14) 

As the structure of a sediment becomes more complex and the 
current paths become more devious, the electrical resistivity will 
increase. It is at this point that the knowledge of both the elec- 
trical and hydraulic properties of sediments and the relationship be- 
tween them becomes important if physical properties are to be proper- 
ly interpreted by electrical resistivity measurements. Much work has 
been done on the flow problem through porous media. The general flow 
formula, be the flow fluid, electrical, thermal or chemical, is: 



F = L* X« A 



(15) 
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F is the flow rate, 

X is the gradient, 

A is the cross-sectional area, 

L is a coefficient of conduction. 

Numerous equations directly inter-relating hydraulic, thermal, and 
electrical properties of porous, solution-saturated materials have 
been derived. For example, Mutt and Berg^ developed prediction 
equations by investigating thermal and electrical conductivities of 
sandstone rock and ocean sediment and relating the two conductivi- 
ties through the common parameter of porosity. Although such a sys- 
tem seems reasonable it is important to remember that the success of 
any electrical method for measuring thermal conductivity in a sedi- 
ment would depend on the range of electrical conductivity of the 
soil particle encountered in a given situation. Therefore it is 
imperative when inter-relating flow problem that these differences 
be understood. 

Schopper^ contends that if quantitative evaluations between 
electrical and hydraulic properties are to be improved, the differ- 
ences between electrical and hydraulic tortuosity should be taken 
into consideration. The essential difference is that hydraulic 
tortuosity takes into account a direction factor in addition to the 
normal tortuous path length, while the normal tortuous path length 
is the only factor involved in electrical tortuosity. Hydraulic 
permeability is a good indicator of hydraulic tortuosity, and there- 
fore improved relationships between porosity and Formation Factor 
would take into account hydraulic permeability. 
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Fatt and Wyble and Glanville investigated the effects 

of applied pressure on electrical resistivity values of reservoir 
rocks and sandstone. These rocks of low porosity were placed under 
very high pressures ranging from 1,000 to 12,000 p.s.i.g. The ef- 
fect of pressure on electrical resistivity seems to be much greater 
than the effect on porosity, since electrical resistivity was 
increased much more than porosity decreased. Similar work on sedi- 
ments appears to be nonexistent, 

Ma thematical Models 

When selecting mathematical models which best fit a porous 
medium such as a saturated sediment, it is necessary to be aware of 
the influence of the charge on clay particles, clay type, salinity 
of the continuous phase, sediment particle size, shape, distribution, 
and structure. In order to further refine models, awareness of other 
factors such as hydraulic permeability and consolidation pressure 
may be useful. In addition the assumptions of theoretical models 
need special consideration, and extrapolation beyond the range of 
experimentation (in empirical models) is not recommended. 

The electrical conductivity of dispersions has been studied 
extensively. Some studies date back over a hundred years. Only 
those theories that are most likely to apply to sediments are con- 
sidered below, 

21 

In 1881 Maxwell conducted one of the earliest theoretical 
investigations and derived an equation for the assemblage of spheres 
in a continuous media. The assemblage was a suspension of spheres 



that were not in physical contact with each other. Maxwell's 
equation, in non-dimensional conductivity terms, follows: 
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K 

m 



3 K d - 2 nK d + 2 n 

K.n - n + 3 
d 



(16) 



If the spheres are assumed to be non-conductive , i.e. Kj approaches 
zero, Maxwell's equation reduces to: 



K 



m 



2 n 

3 - n » 



(17) 



or in resistivity terms: 



FF 



R 



m 



3 - n 
2 n * 



(18) 



This equation can be applied to a sediment mixture consisting of a 
fluid of relatively low electrical resistivity and a network of soil 
grains, not in contact with each other, that are insulators of high 
electrical resistivity. For instance soil granules composed of the 
following common minerals shown in Table 2 would have a high resis- 
tivity. 



TABLE 2 



ELECTRICAL RESISTIVITY OF SOME COMMON SOIL MINERALS 



Mineral 



Resistivity (ohm meter) 



Clay minerals 
Quartz 



10 U - io 12 
io 12 - io 14 



io 7 - io 12 



Calcite 
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When these soil granules of high resistivity (r^) are placed in a 
fluid such as sea water which has a relatively low resistivity (r ) 
(0,25 ohm meters), the non-dimensional resistivity of the dispersed 
phase (soil granule) can be expressed as: 




> CO, 



(19) 



or the inverse relationship: 

k H 

K d = IT > °- (20) 

c 
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Wiener and Woodside and Messmer developed the Series and 
Parallel Model Equations which are considered the limiting relation- 
ships for porosity determination by means of conductivity, A 
complete development of these equations is given in Appendix A, 

In the Parallel Equation, the interstitial water and solid 
soil particle are assumed to be parallel to current flow. The re- 
sulting equation in non-dimensional form follows: 

(Parallel) K = K, + n(l - K,). (21) 

m d d 



The Series Equation evolves when the solid soil particles 
and interstitial water are assumed to be perpendicular to the flow 
of current: 

K d 

(Series) K = - — ; rr. 

'm 1 + n(K, - 1) 
d 



( 22 ) 
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Wiener recognized that other models would fall between the Series 

and Parallel Equations. He proposed a linear combination of the two 

models. Later this linear combination was referred to as the 

23 

Three-Resistor Model by Woodside and Messmer. The Three-Resistor 
Model is given below: 

= p (Parallel Model) + (1 - p) (Series Model). (23) 

In this combination p represents the fraction of the sediment 

described by the Parallel Model. 
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Woodside and Messmer also proposed a Geometric Equation 
in which the sediment is assumed to be the geometric mean of the 
conductivities of the components of the sediment: 

K = K, (1 " n) . (24) 

m a 



However if it is assumed that the solid soil particle is non- 
conductive, then approaches zero. Consequently the Series 
Equation and the Geometric Equation have no meaning, and the 
Parallel Equation is reduced to: 



K 

m 



= n > 



(25) 



and the Three-Resistor Model is reduced to: 

K » pn. 
m K 



( 26 ) 



An extensive review of the work up until 1959 including 
investigations by Rayleigh and Runge was conducted by Meredith.^ 
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He concluded that effects of particle size, shape, distribution and 
structure cannot be ignored if accurate theoretical predictions of 
dispersions are to be made. Meredith developed the Distribution 
Model for treating a concentrated dispersion of spherical and non- 
spherical particles. The relationships in the Distribution Model 
that are most applicable to a sediment environment follow: 



nonconducting 
spheres : 



, _ 8n (n + 1) 

m ~ (5 - n) (3 + n) 9 



(27) 



nonconducting 
randomly oriented 
rods : 



, _ 9n (n + 1) 

m (4-n)(5+n)* 



(28) 



rods perpendicularly 
oriented to field 
current : 



K 



m 



n (n + 1) 
3 - n 



(29) 



Although the above relationships are presented in non-dimensional 
conductivity terms, it is simple to convert them to Formation Factor 
by using the relationship: 



FF 




(7) 



Formation Factor-porosity relationships proposed by Maxwell, Wiener, 

Meredith, Woodside and Messmer are presented in Figure 3. Wherever 

K, is used in formulation it is assumed to be 0.01. Since this value 
a 

is assumed to be very small, 0.01 was arbitrarily chosen. 

Electrical resistivity methods became established as a Geo- 



FORMATION FACTOR 
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POROSITY (%) 



Figure 3* Formation Factor-Porosity Relationships Proposed by 
Maxwell, Wiener, Meredith, Woodside and Messmer. 

(In the Series and Geometric equations the value of 
~ 0.01 is assumed; in all other equations - 0. ) 



S3IU3 
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physical exploration technique in 1927 when Conrad Schumberger 
showed that electrical resistivity measurements of an oil well 
(then called electrical coring), could distinguish between produc- 
tive formation rocks and those that are non-productive. This 
development stimulated extensive research on electrical resistivity 

measurements as an indicator of physical and chemical properties. 
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In 1941 Archie defined formation resistivity factor i.e., 
Formation Factor, as: 



FF = , (30) 

c 

where r is the resistivity of the saturated sediment and r is the 
m J c 

resistivity of interstitial water. Archie also empirically deter- 
mined an equation relating FF to the formation porosity: 

FF = n“ m , (31) 

where n is the porosity of the formation, and m is the slope of the 
line in a log-log plot of Formation Factor versus porosity diagram. 
From his experimental work, Archie determined m to be in the range 
of 1.8 to 2.0 for consolidated sandstones, and about 1.3 for loosely 
consolidated sands • 

Winsauer et al.^ investigated FF to n plots and suggested the 
relation should be: 



-m 

t 



FF =* an 



(32) 
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where (a) is not always unity, but is often a factor accounting for 
the amount of cementation between mineral grains. Therefore (a) is 
referred to as the cementation factor. From equation (32) the 
"Humble Formula" was developed and is probably the most widely used 
equation in the well logging industry: 

FF = 0.62n" 2,15 . (33) 

A comprehensive review of geophysical well logging through 1961 has 
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been compiled by Dakhnov. He summarized the factors that affect 
electrical resistivity of a porous media by a general formula. If 
the equation is applied to a saturated sediment, it becomes: 

r m = f l (c) f 2 (n) f 3 (S) V T) f 5 (Q) f 6 (r d ) W’ (34) 

where r is the resistivity of the sediment, 

c the amount of clay/silt in the sediment, 
n the porosity of the sediment, 

S the partial saturation of the sediment, 

T the temperature of the sediment, 

Q the cation exchange capacity of the soil mineral, 

r^ the resistivity of the soil mineral, 

r the resistivity of the interstitial water, 
c 

Since some of the variables are sometimes omitted, the complexity of 
the above equation is often reduced. For example if the soil miner- 
al (r^) is nonconduc tive , this term is omitted. In a saturated 
sediment (S) is a constant equal to one. If the environment is 
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saline the interstitial water reduces the effect of cation exchange 
capacity (Q) and cation exchange capacity is usually ignored. How- 
ever if the interstitial water has a low ionic concentration this 
term must be considered. The temperature (T) can usually be measured 
and a correction to electrical resistivity can be made based on the 
following formula: 

r 25 - r T [1 + t c (T - 25°)], (35) 

where is the resistivity at 25° Centigrade, 

r T the resistivity measured at temperature T, 
t the temperature coefficient of resistivity, which usually 
is about 2.5 percent per degree Centigrade. 

In most situations then, the resistivity of a saturated sediment 
(at a set temperature) will be determined by three parameters if the 
environment is saline: r^, c, and n; and by four parameters if the 
ionic solution of the interstitial water is low: r^, c, Q and n. 

Dakhnov also obtained a theoretical relation between Formation 
Factor and porosity of uniform nonconducting spheres in a conducting 
fluid: 



FF 



1 + 0.25(1 - n) 



1/3 



1 



(1 - n) 



2/3 



(36) 



Semenov suggested a theoretical equation for cemented rocks 
consisting of grains which have the forms of ellipsoids of revolu- 
tion. The formula for the Formation Factor along the axis of 
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revolution is: 



In n 




In (1 - v ) 
o 



(37) 



In this equation f is defined as: 



3 

e 



(38) 




where 



e 




(39) 



e is the eccentricity of an ellipse, a and b being semimajor and 
semiminor axis of the ellipsoid. The extreme cases of this equa- 
tion are when f = 0 (nonconducting grains of a flat disc shape) 
and when f = 2 (nonconducting grain would be spherical). This 
theoretical work, although developed for the electrical resistivity 
of rocks, combined two additional important factors, that of parti- 
cle distribution which is termed the sorting coefficient ( v Q )t and 
a particle shape factor (f). These two factors would also be 
important in a sediment environment. 

Formation Factor, porosity relationships proposed by Archie, 
Dakhnov and Semenov are present in Figure 4, The Humble Formula 
is also presented in Figure 4, 



FORMATION FACTOR 
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Figure 4. Formation Factor-Porosity polat ionships Proposed by 
Archie , Dnkhnov, Semenov and the Humble Formula.* 



32 



More recent works which are not reviewed previously in the 
literature are discussed below. These recent works are either di- 
rectly applicable or actually conducted on sediments, 

Atkins and Smith^ have experimentally found m, which they 
call shape factor exponents, for various mineral particles: 



Mineral Particle 


m 


Na Montmoril loni te 


3.28 


Ca Montmorillonite 


2.70 


Illite 


2.11 


Kaol ini te 


1.87 


Sand 


1.60 



They also concluded that (a) the cementation factor in formula 32 
was not unity but was a factor of proportionality which accounts for 
amount of clay to sand ratio. Only if a system consisted of one 
type of mineral particle would (a) be unity. 

It is interesting to note that if m could be found by experi- 
mentation, then for that type of clay only the determination of the 
Formation Factor would be necessary and the porosity could be esti- 
mated. Although this study was conducted in order to understand 
how the clay mineral fraction influences reservoir rocks, the results 
could be very useful in the prediction of porosity from Formation 

Factors in sediments, 

9 

Boyce took 43 electrical resistivity measurements on sedi- 
ments taken from cores obtained from seven locations in the Bering 
Sea, some of which were 600 nautical miles from each other. Using 
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a consolidome ter he squeezed the interstitial water out of these 
sediments and measured the electrical resistivity of this water in 
order to obtain Formation Factors. The porosity of these measure- 
ments ranged from 58.3 to 87.4 percent. By the method of least 

\ 

squares for an exponential function the following formula was 
derived ; 

FF = 1.30n" 1#45 . (AO) 



Boyce found that this formula had an error of + 15 percent. He 
stated that this rather large error could be caused by different 
types and percentages of electrical conductors, clay minerals, and 
pore textures. The large error suggests that 600 miles is too large 

an area to attempt to define with one formula. 

12 

Kermabon et al. studied the electrical resistivity of deep 
sediments from 21 long cores taken in a 4 by 10 square mile area in 
the Tyrrhenian Sea. The porosity of these sediments ranged from 
50 to 87 percent. By assuming the resistivity of the interstitial 
water was the same as the bottom sea water, Formation Factors were 
determined. More than 2500 measurements of porosity and Formation 
Factors were compared. A third degree polynomial curve was found 
to best fit their data relating porosity and Formation Factor 



n(7.) = -5.9021 FF 3 + A0.0A16 FF 2 - 105.3889 FF + 171.250A , (Al) 



or, 



FF = 



1 + 0,7193 (l-n 1 ^ 615 ) 



n 



1.A615 



(A2 ) 
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According to Kermabon et al. the above equation fit the data 

obtained with only a Hh 2.5 percent porosity prediction error. 

13 

Smith studied electrical resistivity measurements of North 
Atlantic deep-sea cores and in situ resistivity measurements on the 
continental shelf around Wales, England. The resistivity of the 
interstitial water was either estimated from bottom water salinity 
or actual core water salinity measurements were taken in the labora- 
tory. A best fit straight line through 99 measurements resulted in 
the following equation: 

FF = 1.35n“ 1,2 . (A3) 

Smith noted that the data clearly grouped into two classes: for 

silts and clays of porosities greater than 60 percent, 

FF = n" 2 , (44) 

and for sands and coarse silts of porosities less than 60 percent 

FF = n -1,5 . (45) 

He concluded that there was no need to use a polynomial formula of 

12 

the type described by Kermabon et al., but that once a Formation 
Factor is measured an approximate value of porosity could be ob- 
tained from a general empirical equation, and that after selective 
sampling a better value for porosity could be obtained by a more 
refined equation. 

Porosity versus Formation Factor equations are plotted in 
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Figure 5 from Boyce, Kcrmabon et al., and Smith. Figure 6 shows 
Atkins and Smith’s formulations using equation (32) for sodium 
montmorillonite (m = 3.28), illite (m = 2.11), kaolinite (1.87) 

and sand (m = 1.60). 

26 

Pautot was able to correlate several sedimentary layers by 
taking electrical resistivity measurements on sediment cores ob- 
tained some distance apart. He did this work both in a lake and in 
a marine environment. Although no formulation between electrical 
resistivity and sediment properties was given, the ability to match 
horizontal layers may indicate a similarity between physical and 

chemical properties with electrical properties of sediments. 

27 28 

Bouma et al. and Chmelik et al. developed laboratory and 
in situ electrical resistivity equipment for marine sediments. 

Their objective was to make correlations on lithology and geotech- 
nical properties. Although no formulations were given, they at- 
tempted to correlate electrical properties with the chemical, phys- 
ical and engineering characteristics of sediments. Properties such 
as pH, Eh, water content, density, carbonate content, grain size 
analysis, X-ray radiographs, photographs, cone penetrameter and 
vane shear measurements were compared with electrical resistivity 
measurements. Certain trends were apparent from their data. Elec- 
trical resistivity was indirectly proportional to the water content 

and percentage of clay size present. 
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Ball studied the literature in order to determine the 
feasibility of using electrical conductivity techniques for sediment 
porosity prediction. He concluded that the following points need 
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POROSITY (%) 



Figure 5. Formation Factor-Porosity Relationships Propose*! by 
Doyco, Kermabon, et. al. r and Smith. 
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POROSITY (%) 



Figure 6. Formation Factor-Porosity Relationships Proposed by 
Atkins and Smith for sodium Montmorillonite, TlUte, 
Kaolinito and Sand. 
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further work both in the laboratory and in the marine environment* 

1* What theoretical or empirical equation best describes 
the electrical conductivity-porosity relationship? 

2. What is the effect of clay on electrical conductivity 
measurements ? 

3* How can the conductivity of interstitial water in marine 
sediments be conveniently determined? 

4. What is the best method for restricting the volume over 
which a measurement would be made? 

The first two questions are considered in this section of the paper, 
while the last two questions are investigated in the second section* 
The main motivation behind this literature review and the following 
laboratory experiments is to justify the use of electrical resisti- 
vity techniques as an in situ measuring method to predict porosity 
in a marine environment* Consequently sediment type, interstitial 
water salinity, and formulations applicable to the marine sediment 
environment were studied more extensively than those applicable to 



a fresh water environment 



EXPERIMENTAL PROCEDURE 



Materials 

The following types of materials were obtained for testing: 

1. Kaolinite clay (Albion Kaolin from Babcock and Wilcox 
Refractories Division) , 

2. Illite clay (Grundite from A. P. Green Refractories 
Company), 

3. Montmorilloni te clay (Black Hills Bentonite from 
International Minerals and Chemical Corporation), 

4. Providence silt, 

5. Ottawa sand (well-rounded grain shape), 

6. Two types of glacial sands (angular grain shapes), 

7 • Four types of marine sediments. 

One type of marine sediment was obtained by a Ewing gravity 

corer and the Deep Ocean Sediment Probe (DOSP) in April 1970 from 

shallow water near Ponce, Puerto Rico. The properties of this sedi 
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ment were extensively investigated by Lewis. The other three 
types of marine sediments were taken by an Electrical Resistivity 
Measuring System. These samples were taken during the Fall of 1971 
from Narragansett Bay. 

The properties of the marine sediments and the other labora- 
tory materials are presented in Table 3. Co-ordinates of marine 
sediment cores measured in this investigation are given in Table 4. 



